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a ^ IAW~ STABILIZATION OF LANDING CRAFT 
RICHARD JON HINKLE 

Submitted to the Department of Naval Architecture and Marine 
Engineering on May 20, 1966 in partial fulfillment of the 
requirements for the Master of Science degree in Naval Archi- 
tecture and Marine Engineering* 

The purpose of this work was to investigate the applica- 
tion of pump jet propulsion and yaw angle stabilization for 
a landing craft going through the surf* Linear theory was used 
throughout the analysis of the control problem. The stability 
indices and derivatives were determined in part by experiment 
and the remainder by analytic techniques* 

It was found that LCVP landing craft are dynamically 
unstable without a ekeg. Stability could possibly be restored 
by a feedback control system, but the far more practical solu- 
tion is to start with a dynamically stable craft. This will 
provide the craft with a greater margin of stability in all cases* 
Although it is felt that the craft can not be satisfactorily 
us*?d without a skeg, much of the analytic data in Appendix A 
may be used to cover the case with a skeg. Reference (3) con- 
tains the appropriate modifications needed for the changes. 

Thesis Supervisor: MARTIN A. ABKOWITZ 
Title: Professor of Naval Architecture 
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velocity of hull ceaiter of gravity in the forward direction 
differential operator 

velocity of hull C.G. in starboard direction 

longitudinal distance from L.C.G. to the C.G. of the 

lateral added mass 

distance to L.C.G. from the reference point for forces,, 
moments, velocities etc. 

Longitudinal distance from L.C.G. to the center of 
pressure at which the lateral force f acts (taken as 
center of area of hull profile) 

approximated by half the prismatic coefficient 

force in the starboard direction 

mass density of water 

yaw angle 

subscripts meaning partial derivative with respect to 
r&v respectively 

subscript meaning maximum value 



INTRODUCTION 



Landing craft used in amphibious warfare are susceptable 
to broaching while transiting the surf. The crait, in the surf, 
is particularly vulnerable because it is in a following sea 
and also the relative velocity of the water in the vicinity 
of the rudder is small, making the rudder less effective© 

A possible nmthod of reducing this broaching problem 
is to use a yaw angle control system. By replacing the screw' 
propeller and rudder with a pump jet with a controlled dis- 
charge angle, the problem of decreased rudder effectiveness 
would be eliminated. A pump jet propulsion system of the 
sire required for the craft needs some means of power assist 
to vector the jet stream for steering. This vectoring system 
might easily be converted to an automatic steering control 
system to minimize the yaw angle and eliminate broaching. 

Tb^ economics of war require that in order to be effective 
the craft must be cheap so the control system must be of 
relatively simple design and aceuate for small peri r •? .■* of 
tire such as transiting a surf line. 

The equations of motion for the craft are a pair of 
linearized second order differential equations relating yaw 
and sway to the yawing moments and the transverse force. 

The coefficient© for the various terms in these equations 
were determined by either experiment on a model LCVP or the 
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analytical approach of Jacobs (2). The current LCVPs are 
equipped with a skeg for improved dynamic stability and 
protection of the propeller* The tests run on the model at 
the towing tank were all conducted without this skeg* 

The analysis of the control problem was done using 
linear control theory. The criteria of stability for the 
system was either that of Nyquist or by root locus. A model 
of the system for the analog computer was made, but due to the 
instability the results are not shown* 



k 



PkQCEDUKE 



From the general linearized equations of ship motion of 
Abkowitz (1) the equations for yaw and sway were obtained by 
decoupling the equations for the other modes of freedom# The 
resulting equations are: 



( 1 ) 



( 2 ) 



b ll b 12 , v b 13 dv b l4 

(Y -mx r ) — ~ +• (Y -mu )— -+ (Y -m)~ + Yv=-Y - Y , 
r G ,.2 r o , . r dt v rud d 
dt dt 

(N -I + (N -mx^u + (N + N v = -N ,~N, 

* * r ° ° dt » ° dt v rud d 

b 21 b 22 b 23 ‘24 

where Y , and N , stand for the forces or moments caused by 

the steering device and the subscript d for the disturbance# 

Noting that any force can be resolved into a force and a moment 

at some other location, the force, due to deflecting the pump 

jet, will become the Y r ) ( and N ^ at the reference position 

for the equations. ^ rU( j i® simply Y d multiplied by the distance 

k between its place of application and the reference point# A 

similar treatment applies to disturbance forces and moments# 

deducing the pair of simultaneous equations to one equation in y' 

and Y force we haves 



^ B b 2lV b n b 23 )a * lb l4 b 21 tb 13 b 22- b 24 b 14- b 23 b 12 )s2 * 



lb 14 b 22“ b 24 b 12 , 3“ J n«d(2 D 23' kb 13^ S * ^ b 24" kb l4?] * J dB 



V k 'V a < 



(b 24- k ' b 143 



' i ” " _ V C,. „ :-i •ftifflt --c-'-i 
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Thv block diagram representation of the ship i6 shown in figure i« 

The determination of the stability derivatives was done by 
experiment where the necessary equipment was available and the 
remaining derivatives evaluated analytically. The evaluation of 
N v was done by towing the model in the tank at a small yaw angle 
and measuring the moment. The moment was measured by construct- 
ing a heave rod with the lower eight inches machined round and 
instrumenting the round portion with a full bridge of strain 
gauges. The full bridge array allowed the measurement of torque 
without any effect due to the bending moment in the bar. The 
output of this dynamometer was averaged by an integrator to ob- 
tain the result. Measurement of Y f was accomplished at the same 
time by mounting the thrust block athwartships in the model to 
measure the force in the Y direction. Buns were made at three 
separate yaw angles at each of four different speeds. These 
measurements were made at the longitudinal center of gravity in 
order to have several terms drop out of the equations. The 
remaining stability derivatives were calculated using the method 
set forth in Estimation of Stability Derivatives -■ 

Various Ship Forms t and Comparison With Experimental Results by 
Jacobs. This method employs a strip theory technique with 
lateral added mass coefficient along with Lamb’s coefficients 
of accession to inertia and hydro foil, theory. 

The use of pump j-oi propulsion introduces an additional 
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FIGURE I 



Block Diagram of Ship 
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f» .tor in both the Y and N derivatives. Due to the relative!./ 

r r 

high velocity of the water jet, a rotational momentum change force 
become significant. The force derivative is equal to twice the mass 
of the water in the duct times the velocity of the water jet. The 
moment derivative is derived by applying this force half way down 
the longitudinal length of the discharge duct. For purposes of 
calculations a ten foot duct with a six inch diameter was assumed. 
The thrust of the pump jet was determined by scaling up model 
resistance data obtained in the towing tank, an atter.pt v&l 
made to gather resistance data for the model in a surf but the 
waves would not break for a long enough distance to gather 
any data. 

Because pump jet propulsion was employed, the skeg was not 
needed for protection of the propeller and was not on the model 
during the tests. The elimination of the skeg made the model 
dynamically unstable in straight line motion. 

The automatic control system was chosen to be one of few 

parts* The quantities f and are measured by a free gyro and a 

r* gyro. It may be acceptable to have accura/r 

periods of time such as when going through the surf. In this case 

the measurement of f can be accomplished by using the gyro as a 

free gyroscope rather than as a null device. If the gyroscope is 

fioe, some means of sensing its position without applying much 

tcrque, such as a differential transformer, will have t® be used to 

« 

avoid unwanted precession. The values of & and f times their appro- 
priate gains are summed and compared with the rudder angle and 
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the difference applied to a push pull amplifier. The output of 

the amplifier is a solenoid shown in figure II. By keeping the 
magnetic field in the solenoid large enough the plunger will 
become magnetically saturated. The force exerted by the coils 
on the plunger will then become directly proportional to the 
current. The position of the plunger is used to control an under- 
lapped spool valve of the constant pressure hydraulic system. 

The reaction of the spool valve as a spring and damper were con- 
sidered so that this portion of the system could be critically 
damped. This was done to help reduce the high fr»-r»ency oscilla- 
tions of the vectoring duct. The output from this four way valve 
drives the rams that position the jet vectoring duct. 

The removal of the skeg made the craft dynamically unstable. 

The next step was to investigate what could be done by the con- 
trol system to restore the dynamical stability. First, the effects 
of control, proportional only to yaw angle, was investigated. 

It was shown by the Nyquist Stability Criteria that, at the speed 
chosen for investigation, no amount of gain would suffice for 
stabilization. Proceeding to a control proportional to yaw angle 

\ 

and yaw rate again, it could be shown that the system would be 
ineffective. Using a root locus plot, it can be seen that the 
right most two poles would continue to stay on the right half plane 
for all values of gain. By studing the root locus further, it 
can be seen that a system proportional to yaw plus the first and 
second derivatives would be the minimum requirement to produce 
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a dynamically stable craft by adding two complex c n ,,ga* zer ~ 
to the root locus through the use of appropriate feed back gains, 
the unstable roots could be brought over to the left half plane. 
The design of such a system was dismissed however, kydill (3) 
claims that such systems are not satisfactory due to excessive 
high frequency oscillations of the steering device. Also, the 
expense of the extra equipment needed for the control isould be 
eliminated by reinstalling the skeg and making the craft dynam- 
ically stable to begin with. 
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The results of this investigation were based on the following 
assumptions regarding the design of the system; 

a. longitudinal length of jet duct is 10 feet 

b. jet duct diameter is 6 inches 

c« Y force from deflected jet acts 15«9^ feet from L.G.G. 
d. density of the hydraulic fluid 57 lbs/ft^ 

2 

3 

e« rnaoo of spool valve and plunger is 6*88 x 10 ^ 

to discharge coefficients for spool valv*» ' 
g» width of spool valve ports are 0«?5 inches 
h. the rams that drive the jet vectoring conduit are 1 
square inch in area ana are located 6 inches from the 
swivel, therefore the discharge angle is equal to the 
ram travel in inches over 6 
i* supply pressure is constant at 1000 psi 
j. the radius of gyration equals one forth the length 
The thrust to be developed must equal the resistance of the 
or Aft for constant velocity operation* The model was tested in the 
towing tank and the model resistance data and scaled up result 
are tabulated in figure III. A plot of the full scale resistance 
is shown in figure IV. 

The qualities Y and N were found experimentally in the 
towing tank. The three yaw angles used during the testing were 
3.12°, 5*92°, and 9 » 30 ° and the model speeds varied from 1.3 knots 
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to 2.9 knots. Tills corresponded to speeds between *4.75 and 10.5 
knots for the full scale craft. A table of the experimental data 



is presented in figure V. The data was non-dementionalized by 

Q 2 2 P 2. 3 

dividing forces by r / u 1 and moments by y u 1 . Both N and Y 

X v v 

were found to be dependent on Froude Humber as is shown by the 
plot in figure VI. The variation of yaw angle due to the twist- 
ing of the rod was computed for the largest torque experienced 
and the effect wj*s found to be insignificant. 

The data for the remainder of the stability derivatives and 
indices was computed by the method outlined in reference 2, 

The method of non r dementionalizing used here was to divide forces 
by U / r The length used in all calculations here was that of the 
waterline 3^*36? feet. The actual computations of these deriva- 
tives and indices are contained in appendix A. 

The value of stability indices, derivatives, and resistance 
for speeds of ^ and 7«5 knots are listed in figure Vll. 

The spool valve is an underlapped four way valves with 
ports 0.75 inches wide and discharge coefficients of 0.7 at each 
port. The equations for flow rate for such a valve is that 



2k 



where j is ram rate, X is spool value displacement. 



P is supply pressure, k is the product of discharge coefficient, 
s 

port width and the square root of two times the acceleration of 
gravity divided by the mass of the density of the hydraulic flui<i 0 
This equation gives j=3&0Qx for this system. The 
spring reaction force of the valve is given by the equqtions 
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k=2 x discharge coefficient x port width x x cos 69° 
or k=2 x 0.7 x 0.75 x 10 3 x 0.358 = 376 1 ? 



The damping reaction of the spool valve is expressed by the 
equation; 

C=discharge coefficient x port width 

where L is the difference between distance from supply port to 

load port and load port to drain port. With L equal to 2, the 

value of C is 0.7 x 0.75 x 0.293 x 2=:0.308 — 

lbs 

By using a spring with a constant of 24 __ to center the 

• lb* 3 

plunger, the total spring constant becomes 400 v~ --» The plunger~ 
spool valve combination is assumed to have a mass of 6.88 x 10 3 ^ -^T, s 5 ' 



m 



To critically damp this system, a damping constant of 3»32 
is needed. An additional damper with a constant of 3«0T 
would be needed to critically damp the system. 



lbs~see 



m 
lbs~sec 
Tn 



The combined control system and ship block diagram for 
speeds of 4 and 7*5 knots is illustrated in figures VIII. 
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FIGURE III 
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0.608 


0.288 


4.55 x 10 5 


0.05 


0.00974 
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rA 
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f— 1 


4.53 x 10 3 


1.116 


0.528 


7.91 X 10 3 


0.1575 


0.00914 


4.615 
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10 3 


4.52 x 10 3 


1.489 


0.705 


1.06l x 10 4 


0.283 


0.00924 


4.363 


X 


10 3 


4.88 x 10 3 


2.016 


0.954 


1.428 x 10 4 


0.668 


0.01012 


4.119 


X 


10 3 


6.00 x 10 3 


2.467 


1.168 


1.750 x !0 4 


1.293 


0.01533 


3.969 


X 


10 3 


1.116 x 10 2 


2.978 


1.408 


2.115 x 10 4 


2.413 


0.01960 


3.831 


X 


10 3 


1.577 x 10 2 


3.219 


1.522 


2.282 x 10 4 


3.593 


0.02515 


3.779 


X 


10 3 


2.137 x 10 2 





v/j r 


n r 






f 


C Resistance 

1 lbs 
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0.288 
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7.97 x 10 3 
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7.69 x 10 3 


69.1 


4.2 


0.705 


1.97 x 


10 7 


3.035 


X 


10 3 


7.92 x 10 3 


126.4 


5.68 


0.954 


2.67 X 


10 7 


2.913 
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id 3 


8.91 x 16 3 


260.5 


6.98 


1.168 


3.275 J 


c 10' 


'2.837 
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i.4oo x 16 2 


616 


8.40 


1.408 


3.94 X 


10 7 


2.771 


X 


10 3 


1.854 X 10 2 


1182 


9.09 


1.522 


4.26 X 


10 7 


2.743 
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10 3 


2.411 x 10 2 
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FIGURE V 



Experimental Values of Y and N 

v v 



y=3.i2° 



fps 


‘ N 
lbs 


in lbs 


Y« 


N« 


2.27 


0.0613 


2.3 


6.37 x 10* 


4.32 x 10 4 


2.27 


0.0580 


2.03 


5.88 x 10* 


3.81 x 10 4 


3.03 


0.0833 


3.48 


4.72 x 10* 


3.67 x 10 4 


4.17 


0.1868 


5.15 


5.43 x 10 4 


2.875 x 10 J 






^*5.92' 






2.27 


0.1162 


3.70 


1.173 x 10 3 


6.95 x 10 4 


3.01 


0.1186 


5.18 


6.80 x 10 4 


5.53 x 10* 


3.03 


0.1035 


4.75 


5.86 x K5 4 


5.00 x 10 4 


4.17 


0.4166 


8.40 


1.25 x 10 3 


4,68 x 10 4 


4.54 


0,5307 


10.40 


1.341 x 10 3 


4,89 x 10 4 






^ 9 . 30 ° 




2.055 


0.1414 


2.73 


1.749 x 10 3 


6.26 x 10 4 


2.17 


0.1818 


3.75 


I .83 x 10 3 


7.04 x 10 4 


3.03 


0.3404 


6.30 


1.93 x IS 3 


6.64 x 10 4 


4.17 


0.7405 


12.10 


2.22 x 10 3 


6.74 x 10 4 


5.03 


1.2020 


20.40 


2.4? x IQ" 3 


7.80 x 10 4 
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FIGURE VII 



STABILITIES INDICES AND DERIVATIVES 
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-393.5 


Y 
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Y 
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- 210.5 


N* 
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N. 
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V 


N 
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- 670.0 


- 1 , 258,0 


N 

V 


-2855.0 


- 2 , 765.0 


I 


4$,?00.0 


46,700.0 


m 


2,765.0 


2,765.0 


Coriolua Force 


275*5 r 


1 , 602.0 


Coriolus Moment 


- 2 , 880.0 jp 


- 16 , 750 . 0 , 


Resistance 


115.0 


800.0 
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FIGURE VIII 



Block Diagram of Ship and Control System 
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DISCUSSION- OF KESULTS 



The instability of the system with control proportional only 
to yaw angle is shown in figure IX , a Nyquist Stability Criteria 
plot at 7»5 knots. There is one pole of the open loop transfer 
function in the right half plane. For u nity gain in the control 
system, there are no encirclements of the minus one point, 
therefore there will be one root in the right half plane of the 
closed loop transfer function. increasing the gain, the small 
loop near the origin (shown in the inset of figure IX) will 
expand and eventually encircle the minus o^ie point once. This 
indicates that increasing the gain will only add to the instabil=> 
ity by putting another pole in the right half plane. The Nyquist 
plot at a speed of k knots is not significantly different to 
change the results regarding stability. 

To add to the sophistication of the control system by include 
yaw rate feedback in the system, we can most easily visualize the 
stability problem by a root locus plot. The determination of the 
exact breakaway point requires the solution of a sixth degree 
equation, but the trend of the system can be seen without the 
exact solution. The poles and zeros of the open loop transfer 
function and the asymptotes for speed of k and 7.5 knots are 
plotted in figures X and XI. The effect of proportional plus 
first derivative control is to add another zero to the plot and 
its position would be governed by the ratio of the gains of the 
proportional t© the first derivative terms. With the total 
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system gain at zero, we will already have one pole of the closed 
loop transfer in the right half plane. As the gain is increased, 
the rightmost two poles will come together, merge, and split, 
indicating complex conjugate roots. The roots would follow 
asymptotes originating from approximately the =95*5 point. At 
no time would all the roots be in the left half plane. The placing 
of a zero at the origin by this method of feedback control, could 
only be accomplished with infinite gain on the derivative feed* 
back and so is dismissed as a solution. The only alternative will 
be the speeding up of the control system, such a* by reducing 
the mass of the plunger and spool valve. This would in effect 
move the first pole on the left of the origin further left. 

After the pole passes to the left of the zero, it may become 
possible, with sufficiently high gains, to bring all the roots 
to the left half plane and stabilize the system. 

Proceeding to the next step would be to make the control 
system proportional plus first and second derivative control. 

Here again, the simplest way to visualize the stability problem 
is with a root locus plot. The effect of such a svstf-” would be 
to add two zeros to the plot. By selection of appi-opnate gain 
radios, these zeros will cause the poles on the right half plane 
to cross the imaginary axis at a relatively low value of gain, 
it asymptotes in this case will be at 60° above and below the 
ra 1 axis and on the negative rea2 axis. The point of inter* 
se'.tion of the asymptotes and the real axis will be governed by the 
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real pai 01 position of the two new zeroso The 60° slopes 
obviously indicates that excessive gain will again place the roots 
in the right half plane and unstabilize the system* As was 
mentioned before v systems with second derivative feedback have 
not proved satisfactory for automatic steering control aboard 
ship* The control system tends to produce relatively high 
frequency oscillations of the steering device* A continual ra©ve~ 
went of the steering control would be necessary to stabilize a 
dyni^r all) mutable ship, but these systems tend to cause much 
more •notion than is need *1 f rt r effective control* The introduce 
lion of a *e< omd degree lag network into the control system is 
/*« • 'i i I - irydi i l ( *>) t« »■ 1 reduce these oscilla.tio.os* 
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CONCLUSION 



At- is evident that any attempt at automatic control for yaw 
s' to libation of an LGVF without a skeg, will be difficult and 
uneooriotf ical . Uy using the skeg to restore dynamic stability 
to the craft before the control system is added, a great deal of 
>i' ■ will foe eliminated. The dynamically stable craft will 
no' iraliy ht %a«ier to handle when steering manually. The use of 

■ * to be a practical method of propelling 
in i font i 1 1 mg landing craft and is adaptable to automatic 
( ontrol. 
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RECOMMENDATIONS 



The mosr important recoimnendation can only be to use a skeg 
to make the Lf.V* dynamically stable to begin with* If this 
scheme of testing is to be used again for this model, Jacob's 
method may be used to include the use of the skeg too* 

The testing in the towing tank should include some means 
f obtaining resistance data while in a surf* This will govern 
t tter«ist avjiiahl* for cntrol. I suggest that some ramp be 
\ *fcal i in the tank near the wave maker to generate breakers, 
fnis ramp * nild have the effect as a reef in causing the 

t i oi ome breakers. I feel any attempts to gather data 
* k og the tanks beach will prove futile, 
d a * he short length of surf. 
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APPENDIX 



APiENDIX A 



The following a brief description of the Jacobs’ (2) method 

of determining stability indices* The first strep wss to deter* 

mine the sectional inertia coefficients C by entering figure XVII 

with the values of sectional area divided by local beam times 

locti! dr It 1**' ) and ioc*i b*u& divided by local draft C^T)* 
oh h 

Next, Iamb’s foef t leionfcs of accession to inertia (k^k^k 0 ) 

are found f r«»a figure XViii by entering with The prime in 

the foil wing *qu*tiv&s will denote non- demen t ionali zed quantities* 

tnote all quantities htre are non°dementionalised by the quanti« 

ties not^d in trv chapter on results) The limits on the 

t.* r b*j* to stern* 

* c» *. 



rV c 

fi\ s ^k.. /2 / G h dx. 
£ ^ '>Tm s 



( 4 ) 






a 



( 5 ) 




( 6 ) 



-S 



fc # , 

;I G 2 



y G £. 

”.5 & 



( 7 ) 

( 8 ) 



iiall i.r, 1 1 u«atic coefficieat 



( 9 ) 



.51 



( 10 ) 



N. =k°^!/c h 2 x 2 dx 
r 2 va s 



Y»=m..=k.,m 
r 1 1 o 


(11) 


Y» ==m_x 
r 2 


(12) 


N» =h® Jx 
r 2 


(13) 


K^'o * "i 


(14) 


»r<4 uhl 


(15) 


- X _ 

N° »!<”■) LV 

r 2 11 i 


(16) 


N S =N 9 £uhl 2 
r r'2 


(17) 



Figures XII through XVI are the calculations required to 
derive these stability derivatives and indices for the LCVP 
tested. The integrations were carried out by Simpson 0 s ( rule 
and all areas found by use of a planimeter. 
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FIGURE XII 
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FIGURE XII (continued) 




h 2 xdx=0.935 x 0.995 x 3.14159 x 38.17=111. 



5 



N« =m_x=111.5 
y 2 
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FIGURE XIII 
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FIGURE XIV 
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FIGURE XV 
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FIGUKE XVI 
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FIGURE VIII 




SECTIONAL INERTIA COEFFICIENTS C $ AS FUNCTIONS OF THE 

LOCAL BEAM -DRAFT RATIO b/h AND SECTION AREA COEFFICIENT, 
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